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Abstract 



We discuss the quantization of the energy levels of two-particle scattering states 
in a finite volume in the case of Bloch-type boundary conditions. A generaliza- 
tion of the Liischer quantization condition is obtained that can be used in order to 
calculate the scattering phases, resulting for example by the strong elastic interac- 
l/^ . tion of two pions, at fixed physical momentum transfers on a sequence of volumes 

of growing sizes. We also give a generalization of the Lellouch-Liischer formula to 
' be used to extract the physical K i— > vrvr decay rate below the inelastic threshold 

Tij- ■ from finite volume calculations. The formula is valid up to corrections exponentially 

\ vanishing in the volume. By using this formula the calculation can in principle be 

I performed on different finite volumes of growing sizes in order to keep under control 

the corrections. 
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1 Introduction 

X 

d ' The calculation of the weak matrix elements associated to the non-leptonic 

decay K ^ mi requires a non-perturbative level of accuracy, due to the 
strong nature of the final state re-interactions. Lattice gauge theory is the 
best candidate among the other methods to achieve a phenomenologically 
relevant result for this process but there are theoretical problems that have to 
be addressed before such kind of calculations can be actually performed. The 
obstacles that have to be removed are mainly due to two different reasons. 

The so called "ultraviolet problem" , concerning the construction of finite 
matrix elements of properly renormalized lattice operators, has been addressed 
by many authors and we will not discuss it further in the rest of this paper 
(see [1] for a recent review of this subject). Here we will focus our attention 
on the so called "infrared problem" . 

The infrared problem, formalized in the well known Maiani-Testa "no- 
go" theorem [2,3], concerns the impossibility of connecting the large time 



correlation functions calcTilatcd on a large euclidean volume with the physical 
decay amplitudes. A solution to the infrared problem has been found in ref. [4] 
where Lellouch and Liischer (in the following referred as LL) derived a relation 
connecting the finite volume euclidean K i— > tttt matrix elements with the 
physical ones, up to exponentially small finite volume corrections. 

The key point in the derivation of the LL formula is that the Maiani- 
Testa theorem applies in the case of very large volumes, where the spectrum 
of the two-particle final state is continuous, but ceases to be valid in the case 
of finite volumes where the two-particle energy is quantized. In the last case 
a kaon at rest cannot decay into two pions unless the physical extension of 
the finite volume is such that one of the two-particle energy levels equals the 
mass of the decaying particle. A simple phenomenological analysis based on 
the results of refs. [5,6,7,8] shows that, when the two pions satisfy periodic 
boundary conditions, the volume on which the energy of the first excited state 
of definite isospin coincides with the kaon mass is of the order of 5.5 fm. Such 
a large value of the volume makes the calculation unfeasible from the practical 
point of view due to the limiting present computer power. Furthermore the 
calculation is performed on a single finite volume and the corrections, although 
exponentially vanishing, cannot be quantified. 

The authors of ref. [9] have shown that a set of LL formulas can be derived 
also when the energy of the outgoing two-pion state does not coincide with 
the kaon mass. In this case the resulting infinite volume decay rates are not 
calculated at the physical point, i.e. at the values of the two-pion energies 
fixed by relativistic kinematics. Nevertheless one can repeat the calculation 
for different finite volumes and then extrapolate the results to the physical 
point. In ref. [9] it has been also pointed out that the finite volume corrections 
to the LL formulas could result to be sizable if only few two-particle states 
have energies below the inelastic threshold. 

In this paper we study the quantization condition of the energy of a two- 
particle state on a finite volume by using a generalized set of boundary con- 
ditions (^-BC) for the two particles. In ref. [10] we have shown that the use 
of 6'-BC in framework of lattice calculations makes possible a continuous mo- 
mentum transfer between one-particle states. In the following we obtain a 
quantization condition for the energy of a two-particle state on a finite vol- 
ume that generahze a result previously obtained by Liischer [7] in the case of 
periodic boundary conditions. We also give a generalization of the LL formula 
and show that by using the ^?-dependence of this formula and of our quantiza- 
tion condition it is possible to find a sequence of finite volumes of growing sizes 
on which the calculation of the X i— > tttt decay rates can be performed at the 
physical point. Using the values of the scattering phases predicted by chiral 
perturbation theory we show that the physical amplitudes can be calculated 
on volumes of the order of about 3 fm and, at the same time, the residual 
functional dependence of the results upon the volume can be quantified. 
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The plan of the paper is as follows. In section 2 we set up the notations. In 
section 3 wc derive the quantization condition while the generalized LL formula 
is given in section 4. In section 5 we study the volume dependence of the two- 
pion spectrum for some particular choices of the ^-angle and in section 6 we 
draw our conclusions. Some technical details needed for the derivation of the 
quantization condition are given in the appendices. 



2 Two— particle states in a finite volume 

The spectrum of a two-particle state on a finite volume in quantum field 
theory has been already studied in great detail in refs. [5,6,7,8] in the case 
in which the two particles satisfy periodic boundary conditions. An energy 
quantization condition has been found by establishing in ref. [5,11] a connec- 
tion between quantum field theory and non-relativistic quantum mechanics. 
Indeed, by assuming that the two particles are spinless bosons of equal mass 
m whose dynamics can be described by a scalar field theory of the 0^-type, 
that the reflection symmetry i— > —(f) is unbroken and that the one-particle 
states are odd under this symmetry, an effective Schrodinger equation can be 
written for the two-particle state. In the center-of-mass reference frame this 
equation reads 

- ^Aijix) + ]- f dx' Ue{x, x) ij{x) = Ei{j{x) (1) 

where the parameter E does not represent the true energy of the system, that 
we call but it is connected to the last through 

E = 2vWTmE (2) 

In eq. (1) the parameter /i represent the reduced mass of the two-particle sys- 
tem while Ue{x,x') is the Fourier transform of the modified Bethe-Salpeter 
kernel UE{k,k') introduced in [5]. The "pseudo-potential" Ue{x,x') depends 
analytically on E in the range —m < E < 3m, is a smooth function of the co- 
ordinates X and f decaying exponentially in each direction and is rotationally 
invariant so that one can pass to the radial effective Schrodinger equation. 

Thanks to these observations, in the following we will carry on the calcu- 
lation of the spectrum of a two-particle state in a finite volume by using a 
purely non-relativistic Hamiltonian that, separating the center of mass motion 
from the internal motion, comes out to depend upon the relative coordinate 
f— X — y only 

H^-^A + V{r), r^\\f\\ (3) 
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where A is the Laplacian operator with respect to r. The potential is assumed 
to be spherically symmetric, a smooth function of its argument and of finite 
range, i.e. 

V(r) =0 for r> R (4) 

We will solve the problem on a finite cubical box of linear extension L in each 
direction greater than the potential radius (L > R) and we will assume that 
the potential is periodic of period L 

V{\\f+nL\\) = V{r) for n e (5) 

We can imagine to start from a given finite-range potential V(r) that describe 
the interactions of the two particles and to build a periodic potential as follows 

V{r)= Y: V(||r + nL||) (6) 

nGZ3 

By construction V{r) satisfies the periodicity condition stated in eq. (5). 

There are two differences between the potential V{r) that we are going 
to study and the pseudo-potential Ue{x,x'). The first one concerns the en- 
ergy dependence of Ue but this does not represent a problem because all 
the results we are going to derive will be obtained at fixed E. The second 
difference between the quantum field system and the non-relativistic one is 
that Ue{x, x') does not vanish if one between x and x' is greater than R but 
has exponentially small corrections. Furthermore in the quantum field system 
there are additional exponentially small finite volume corrections that arise 
from polarization effects. For these reasons the results we are going to derive 
in the non-relativistic theory will be valid also in the relativistic theory up to 
exponentially small corrections. 

The matching with a quantum mechanical system could be avoided at all 
following an approach similar to that developed in ref. [9] . In the following we 
will not follow this strategy because, as will be clear later on in the derivation, 
the quantum mechanical analogy allows us to benefit without additional effort 
of a series of theoretical results obtained in the framework of solid state physics 
and useful also in our case. 

3 Quantization condition: lessons from solid state physics 

In this section we derive a powerful relation connecting the energy eigen- 
values of a two-particle state on a finite volume with the infinite volume 
scattering phases of the two particles. By using the result stated in the previ- 
ous section we will perform the calculation in NRQM being the results valid 
also in QFT up to exponentially vanishing finite volume corrections. 



4 



In order to obtain the energy quantization condition wc have to deal with 
the Schrodinger equation for a particle in a periodic potential, i.e. the same 
equation satisfied by electrons, holes and excitons in a periodic crystal. We 
can thus export useful results well known to the solid state physics community 
since a long time, provided that we interpret the "cell size" of the crystal as 
the physical extension of the finite volume, i.e. L. 



3.1 Block's theorem 



In order to better understand the role of our particular choice of the bound- 
ary conditions we want to recall the well known Bloch's theorem 

theorem 1 The wavefunctions of the "crystal Hamiltonian" can be written 
as the product of a plane wave of wavevector k within the first Brillouin zone, 
times an appropriate periodic function: 

Mrl = e''-' uj:{r-) (7) 



where 



2tt 

Uj^{f+ nL) — Uj^{r) and < ki < — (8) 



Let us observe that we are dealing exactly with a sort of crystal Hamiltonian 
(see eqs. (3), (4) and (5)) and so that the result stated from the Bloch's 
theorem apply straightforwardly also in our case. Furthermore, we observe 
that the wavefunctions ipj^{r) do not satisfy periodic boundary conditions but 
the more general set of b.c. 

^g(f+nL) = e^^^-"^ ^/;g(r) (9) 



These boundary conditions have been recently considered in the framework of 
lattice QCD in ref. [10] under the name of ^-boundary conditions, precisely 
in the form 

^g{r + nL)^e'^-'^ ijgir) (10) 



in order to handle in lattice simulations physical momenta smaller than the 
allowed 

p^-^n, neZ3 (11) 

Ij 

in the case of standard periodic boundary conditions. The matching between 
the formalism of ref. [10] and the condition stating the fundamental result of 
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the Bloch's theorem is obtained by identifying 

k^^ (12) 



3.2 Korringa-Kohn-Rostoker theory 

Another fundamental result that we can gain from solid state physics is the 
computational framework developed independently by Korringa [12], Kohn 
and Rostoker [13] and known as the KKR method or the Green's function 
method. A straightforward application of this method will allow us to derive 
the two-particle state energy quantization condition in a simple and (in our 
opinion) clear way. 

The KKR method can be applied under the hypotheses of a so called 
"muffin thin potential", i.e. a periodic, spherical symmetric potential that 
vanishes after a given distance R for each cell of the crystal. After having 
realized that all these hypotheses are satisfied by our potential defined in 
eq. (6), we can start in reviewing the KKR procedure by considering the 
time-independent form of the Schrodinger equation 

(A + q') i^j;ir) = y'(r)^,-(r) (13) 

where we have defined 

g2 = 2nE (14) 

and we have substituted 2fj,V(r) with V'{r). In order to have a formal solution 
of this non-homogeneous partial differential equation we introduce the free- 
particle Green's function as the solution of 

(^Ar^ + q'^^ g{f - fo; q) = S{f - fo) (15) 
Using the Green's function, the solution of eq. (13) can be written formally as 

^fe(^ = + / g{r- ro;q)V'{fo)^j:{fo) (16) 

— oo 

where (/"^(r) is a solution of the homogeneous equation associated with eq. (13) 
with the additional requirement to satisfy the Bloch's condition of eq. (9). 
In the previous equation the integration variable tq spans the whole three 
dimensional space and not only a period. We will set for the moment ^^(r) = 
and we will come back to the case in which the homogeneous solution is present 
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later on. We rewrite eq. (16) as 

+00 

= / dro g(f- fo; q)V' (ro)i^i;(fo) (17) 



We want to observe that no particular conditions are requested to the Green's 
function in order ipj^{r) to satisfy the Bloch's condition; indeed being ipj^{r) 
present in both the members of the previous equation and being the potential 
periodic, the condition of eq. (9) is self-consistently satisfied. For this reason 
we will not require g{f—fo;q) to satisfy any particular periodicity condition. 
As a consequence we can use for this function the well known result 

1 piq\\f-fo\\ 

^(r-r-'o;?)^--^^ (18) 
471 ||r — roll 

The domain of integration in eq. (17) can be reduced from the entire 
world to a single periodicity cell by introducing the "greenian" of the equation 
defined as 

gnir- fo; q)= Yl e*"^i/(f - fo - nL; q) (19) 

neZ3 

Using the greenian definition together with the Bloch's condition and the 
periodicity of our potential, the formal solution of the Schrodinger equation 
can be rewritten as an integral spanning only a period 

V'fel^ = / dro g^{f-fQ]q)V'{fo)'il)^{rQ) (20) 
period 

but, being the potential identically zero for distances greater than i?, the 
integration domain can be further reduced to a sphere of radius R 

= / dfo gnif- fo; q)V'(fo)Mro) (21) 

Sr 

Now we use the fact that, thank to the Schrodinger equation (13), the previous 
relation can be rewritten in a form suitable for the application of the Green's 
theorem 



'^ki^ = J dfo g^if- fo; q){Af^o + 1 )'^k{^o) (22) 

Sr 

Using the identity 

gAij = ^A^ + V • (^V^ - ijVg) (23) 
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and the Green's theorem we end up with a vanishing surface integral 



/ 



'0 

dSn 



9ki'^-ro;q)^- V'fe(ro) 



dro dro 



= (24) 

ro=R 



The previous equation is the quantization condition we where searching for. 
In the following, by using the expansion in spherical harmonics of both the 

greenian and of the wavefunction, we will recast this condition in a system 
of equations expressing the two-particle scattering phases as functions of the 
energy eigenvalues and vice versa. 

3.3 Partial wave expansion of the wavefunction 

Let us consider the first periodicity cell. For distances greater than the 
potential radius, the wavefunctions V'fe(r) satisfy the free particle equation 
and can thus be written as 

M^) = Y.<'Uk)Ri{r-q)YUro) (25) 

Im 

where cim{k) are coefficients to be determined by using eq. (24) and Yim{0, (j)) 
are the spherical harmonics. The radial part Ri{r,q) can be expressed as 

Ri{r, q) = cos Si{q) ji{qr) - sin 5i{q) ni{qr) r>R (26) 

where Si{q) are the two-particle infinite volume scattering phases, ji{qr) are 
the spherical Bessel functions and ni{qr) are the spherical Neumann functions. 
For later use we report the Wronskian relations satisfied by the Bessel and 
Neumann functions 

s\u5i{q) cos 6i{q) 
[ji, Ri] = 5— R> Ri] = 5— 27) 

where, as usual, the Wronskian of two functions is defined as 

[fix),g{x)] = f{x)g\x) - f\x)g{x) (28) 

3.4 Partial wave expansion of the greenian 

The derivation of the partial wave expansion of the greenian is a rather 

involved mathematical exercise. In order to make our derivation of the quanti- 
zation condition of the two-particle energy as clear as possible we give all the 
technical details in the appendix A and report here below only the resulting 
expression 
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9k{r-fo;q)^qJ2jiiqr) Yi^{r)ni{qro) Y^^ih) 

Im 

+ H jM^) ^im(r) Vim,i'm'{k,q) ji'{qro) Yi^m'ifo) (29) 

Iml'm' 

where r < Tq < nL ^ 0, and we have introduced the structure coefficients 

^lm,l'm'{k, q) = 47r ^ i~''DjK{k, q) CjK;lm,l'm' (30) 

JK 

In the previous definition CjK;im,i'm' are the so called Gaunt coefficients (sim- 
ply related to the Wigner 3-j symbols) defined as follows 

CjK;lm,l'm' = J dn^YjK{k)Yil{k)Yi,^,{k) (31) 

while Dim{k,q) are the so called reduced structure coefficients. The time con- 
suming part of a KKR calculation is given by the numerical evaluation of 
the reduced structure coefficients. For this reason there are many equivalent 
expressions of these quantities some of which are given in appendix B. Here 
below we limit ourself to observe that the reduced coefficients can be written 
as 

/f N dim{0,Q) 

Dim{k,q) = (32) 

where the dim{0, (5)'s are dimensionless quantities (see eqs. (B.28)) and where 
we have used the following definitions 

e^kL and <5 = ^ (33) 

3.5 Generalized Liischer quantization condition 

The generalized Liischer quantization condition, can be now easily derived 
by substituting in the condition stated by eq. (24) the partial wave expansion 
of the wavefunctions given in eq. (25) and that of the greenian given in eq. (29). 
We obtain 

Q^jiiQr) Yi^{r) [ni,Ri] ci^ + 

Im 

+ XI 3i{.(l'r)Yi^{r)Ti^^i,^, [jv , Ri']ci'm' = ^ (34) 

Iml'm' 

that, making use of the Wronskian relations of eq. (27), can be rewritten as 

q cos 6i{q) Cim + Yl ^lm,l'm' siii 6i< {q)ci'm' = (35) 
I'm' 
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The previous one is a linear homogeneous system of equations having as un- 
knowns the coefficients q^. The following compatibility determinantal equa- 
tion 



det 



^im,i'm'{k,q) +q5w 5mm' cot Si{q) =0 (36) 



is the quantization conditions for the energy of the two-particle states on a 
finite volume. This equation can be used to calculate the scattering phases 
once the energy eigenvalues are know (for example form lattice calculations) 
or vice versa to calculate the spectrum of the two-particle state given the 
scattering phases (for example by the analytical knowledge of the interaction 
potential). 

We want to stress again that the results obtained in this section arc valid 
up to exponentially small corrections. In particular we observe that in order 
eq. (36) to be used in quantum field theory the physical size of the finite 
volume has to be large enough to exclude the presence of residual polarization 
effects. These observations apply also to the results of the following sections. 

3.6 Singular solutions 

In our derivation of the quantization conditions we have assumed that the 
solutions of the homogeneous time-independent Shrodinger equation 

(A + q') <f>j^{f) = (37) 
where absent. In order ^^^(r) to be different from zero the condition 



e + 2nn 



(38) 



has to be satisfied for some integer vector n at fixed 6 together with a quantiza- 
tion condition that can be derived along the same lines that we have followed 
to obtain eq. (36). This can happen only on particular volumes L and/or for 
particular interactions. In the following we will not be interested in this phe- 
nomenologically not relevant situation and we refer the interested reader to 
ref. [7] for a detailed discussion of the 9 = case. 

3.7 S-wave interactions 



Prom the phenomenological point of view it is relevant the situation in 
which all the scattering phases can be assumed to vanish except the S-wecve 
one 

Si{q) = , l>0 (39) 
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In this case the quantization conditions simphfy extremely. Indeed we have 
C„o,«.„o=(^/<iOj = ^ (40) 

that substituted in eq. (36) together with eq. (39) gives us 

tan5o(g) = -^=-f-^ (41) 
V47r Doo{k,q) 

or, equivalent ly, 

As we have observed in eq. (32) (and shown in eqs. (B.28)), the right hand 
side of the previous equation is a dimensionless quantity that can be computed 
once and forever; we set 



The definition of 0(^, Q) is completed by requiring the continuity of this func- 
tion by respect the variable Q for each value of 9 and by the condition 

</,(^, g = 0) = , < < 27r (44) 



4 Generalized Lellouch— Liischer formula 



In this section we give a generalization of the LL formula [4] that connects 
the weak matrix element of the decay K — > (tttt)/ (being / the isospin of 
the state) calculated on a finite volume with the corresponding quantity in 
the infinite volume limit up to exponentially vanishing corrections. In the 
following we work under the same hypotheses of ref. [4], i.e. we consider a 
theory of spinless pions and kaons of masses such that the condition 

2m^ < rriK < ^tUt^ (45) 



is satisfied. Wc further assume that the pions scatter purely elastically below 
the threshold for the production of four pions and that the kaon is stable in 
the absence of weak interactions. When the weak interactions are switched-on 
the kaon is allowed to decay into two pions and the transition amplitude is 
given by 

T{K ^ tttt) = A{q) e^^°® (46) 
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where A is real, So is the yS'-wave scattering phase of the outgoing two-pion 
state and q is the pion momentum in the center-of-mass frame 

1 



^,ml- 4ml (47) 



In writing eq. (46) we have assumed the standard relativistic normaUzations 
for the one-particle states together with the LSZ constraints on their phases 



(0| ip^ix) |7r(p)) 



ipx 



(0| Mx) \K{p)) 



Zk e" 



tpx 



(48) 



where (/?7r and (^k ai'c the pion and kaon interpolating fields respectively. 

Let us now consider the same theory on a finite volume of linear extension 
L. The finite volume one-particle states are normalized to unity while their 
phases can be chosen arbitrarily. Let us call 



AL{q) 



w 



(49) 



the transition matrix element on the finite volume with Hw the weak inter- 
actions Hamiltonian. Our generalization of the LL formula is given by 

where 0(6*, Q) has been defined previously in eq. (43). We omit the derivation 
of this powerful relation that follows straightforwardly by repeating the same 
arguments of ref. [4] and is valid under the same hypotheses that led us in 

section 3.7 to derive eq. (42) plus two additional requirements on the outgoing 
two-pion state. It has to be non degenerate and must have the same energy 
of the decaying kaon, in order the resulting infinite volume decay amplitude 
to be computed at the physical point. 

In the LL derivation one has ^ = and the two-pions state of definite 
isospin happens to have an energy equal to the kaon mass only on a certain 
particular volume. In our generalization of their formula the ^-dependence 
of (l){9, Q) can be used in order to obtain a sequence of volumes of growing 
sizes on which the calculations can be actually performed (see the following 
section). Successively, by studying the residual functional dependence of the 
results upon the volume it will also be possible to answer to the questions 
raised in ref. [9] about the size of the finite volume corrections coming from 
the presence of the inelastic threshold. 

In the particular case 9 = the first seven energy levels of the outgoing 
two-pion states are non degenerate, as shown in ref. [7]. Care about the pos- 
sible degeneracies of the outgoing states has to be taken also in the case 6^ 7^ 
for the particular choices of the Bloch's angle used in the calculations. 
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A further generalization of eq. (50) by respect the LL result can be achieved 
following the arguments of ref. [9] where the LL formula is shown to be valid 
for all the states below the inelastic threshold and even outside the physical 
point {q = q). 



4-1 Choice of the boundary conditions in numerical simulations 

Both the quantization conditions of eq. (36) and the generalized LL for- 
mula of eq. (50) have been derived in the reference frame of the two-pions 
center of mass. Although a generalization of the quantization condition to 
other reference frames could be obtained along the same lines of ref. [14], 
when these formulas are used in the context of lattice QCD calculations the 
^-angles corresponding to the up, down and strange quarks have to be fixed 
in such a way that the decaying kaon is at rest and that the total momentum 
of the two-pion states is zero. 

Let us consider, for example, the decay of a neutral kaon in charged pions. 
In order to have the kaon at rest, considering that the anti-quarks have a 9- 
angle opposite to that of the corresponding quarks, one has to chose strange 
and down quarks boundary conditions so that 

9s = 9d (51) 



At the same time, the positive charged pion interpolating field will have an 
overall ^-angle given by 

C+ ^^u-^d (52) 



while for the negative pion one has 6*^- = —9^^+^ so that the total momentum of 
the two-pion system is zero. There is no problem to fulfill these requirements 
in the quenched theory but, when one considers the full theory, some technical 
complications arise. Indeed current simulation algorithms require the product 
of the quark determinants (that depend upon 9) to be non negative and, when 
Wilson fermions are considered, this implies that 

Ou = Od (53) 



i.e. a vanishing ^?-angle for the interpolating fields of the charged pions. Nev- 
ertheless this complication does not arise when Ginsparg- Wilson fermions are 
considered because in this case the determinant of each quark is non negative. 
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Fig. 1. The red points represent —tan(p{9,Q) as function of Q for 9 = 0.0. The 
different branches are the two-particle eigenvalues curves. The straight lines are 
the tangents of the scattering phases of two pions, having isospin / = or / = 2, 
computed in chiral perturbation theory at the kaon mass. 




Fig. 2. The points represent — tan(j){9, Q) as functions of Q for different values of 9. 
The case 9 = 2.5 is not shown to help the eye. The straight lines are the tangents 
of the scattering phases of two pions, having isospin / = or / = 2, computed in 
chiral perturbation theory at the kaon mass. 



5 Volume dependence of a two— pion state energy spectrum 

We are going to study the eigenvalues of a two-pion system on finite 
volumes assuming that the scattering phases for I > 4 are small in the elastic 
region (note that 6i{q) is proportional to at low momenta). Under these 
hypotheses the energy quantization condition takes the simple form of eq. (42). 

In fig. 1 we show the opposite of tan [0(^, Q)] for the particular choice 

6 = 0.0. The points in which this function coincide with the function tandolq) 
at fixed 6 and L represent the two-particle eigenvalues on the given volume 
but, in the following we want to perform a slightly different phenomenological 



14 



analysis: 

• we fix the energy of tlie two particle state so that q — q 

• we compute the corresponding scattering phase So{q) (using the one-loop 
chiral perturbation theory results for these quantities [15,16,17,18,19]) 

• we find the values of Qni^) such that the condition 

t&nSoiq) = -t&n (j){e,Q) (54) 
is satisfied. 

• the volume Ln{0) on which the n two-particle eigenvalue is equal to q is 
computed according to 

Ln{9) = 27r^^ (55) 



Following these steps we are able to find for each value of 6 the finite volumes 
on which a given energy level of a two-pion state equals the kaon mass. In 
ref. [4] the authors have not considered the possibility of having 9 and 
their analysis, that we reproduce in fig. 1, gives the following results 

7 = gi(0)=0.89 Li(0)=5.34fm 

7 = 2 Qi(0) = 1.02 Li(0)=6.09fm (56) 

i.e. the numerical simulations should be performed on volumes of order of 
5.5 fm that, in the case of unquenched calculations, could be too large even 
for the next-generation super computers. In fig. 2 we repeat the same analysis 
for other allowed values of 9. As can be seen a careful choice of the "Bloch 
angle" shifts the position of the first eigenvalue curves at lower values of Q 
and consequently, being q fixed, at lower values of the physical volume. 

The resulting volumes corresponding to our particular choices of 9 are 
given in tab. 1. As can be seen the volumes to be simulated are of the order 
of 3 fm that are accessible on present super-computers. A similar analysis can 
be repeated for all the allowed values of 6* (0 < 6'^ < 2%) in order to obtain a 
longer sequence of finite volumes of growing physical sizes on which the decay 
amplitudes can be computed at the physical point. As already pointed out in 
the previous section, the size of the finite volume corrections to eq. (42) and 
to eq. (50) can be quantified by repeating the calculations on this sequence of 
volumes and by studying the functional dependence of the results upon the 
box extension. 
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e 


/ 


Ql 


Li 


0.0 





0.890 


5.34 


0.0 


2 


1.015 


6.09 


2.5 





0.295 


1.78 


2.5 


2 


0.415 


2.49 


3.0 





0.400 


2.40 


3.0 


2 


0.490 


2.94 


4.0 





0.569 


3.41 


4.0 


2 


0.645 


3.86 



Table 1 

Volumes on which the energy of the first eigenvalue two-pion state equals the 
if-meson mass. / is the isospin of the state. 



6 Conclusions 



In this work we have studied the spectrum of a two particle state on a 
finite volume with Bloch's boundary conditions. 

We have found a quantization condition that generalizes a result previ- 
ously obtained by Liischer in the case of periodic boundary conditions and 
that relates the energy eigenvalues of a two-particle system on a finite volume 
with their infinite volume scattering phases. The formula is valid up to ex- 
ponentially small finite volume corrections and allows the calculations of the 
scattering phases at small physical momenta by calculating, for example on 
the lattice, the energy spectrum of the two-particle state on volumes of the 
order of 3 fm. 

From our quantization condition it follows straightforwardly a generaliza- 
tion of the Lellouch-Liischer formula that connects the finite volume ampli- 
tudes for the decays of a kaon into two pions with the corresponding quantities 
in the infinite volume. We have argued that the decay amplitudes can be ob- 
tained simulating on the lattice finite volumes of different physical extensions 
in order to quantify the size of the residual finite volume corrections and thus 
to obtain results with a phenomenologically relevant accuracy. 
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Appendices 



A Derivation of the partial wave expansion of the greenian 

In this appendix we derive the partial wave expansion of the greenian of 

eq. (29) that we have used to derive the generalized Liischer's quantization 
conditions. We start by recalling the well known Neumann's expansion of the 
Green's function 



^TTllr-f'll 

The previous relation is valid provided that r < r' otherwise one has to ex- 
change r and r' in the right-hand side. We insert the Neumann expansion in 
the expression for the greenian that we report below for clarity 

1 piq\\r-ro-nL\\ 

For convenience we separate out the term with n — Q from the remaining 
terms and we use the Neumann expansion identifying f* with fo in the first 
case and with nL in the second case. We end up with 

1 cos ||r — fo 



A II- -II +Y.Dimi.Kq) jM\\r-fo\\)YUf-fo) (A.3) 
where the so called reduced structure coefficients are given by 
Dim {k,q) ^ 
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= ^ e'^-^^Yim(n)[ni(qnL)-iji(qnL)] -i-^dioSmo (A.4) 

neZ3-{0} 

There are many different, aitfiougfi equivaient from tlie matliematical point 
of view, ways to express the reduced structure coefficients some of which are 
more convenient than the previous relation for a numerical computation of 
these quantities. In the appendix B we derive an expression suitable for the 
numerical calculation. 

We need to recall another well known identity that can be easily proved 
by using the expansion in spherical harmonics of a plane wave (see eq. (B.7)) 

i''''jj{qK)YjK{R)^ 

= 47r ^ i^~^'CjK;im,i'm' ji{qr) ji'iqro) Yi^{f)Yi,^,{ro) (A.5) 

lm,l'm' 

where we have called R = ||r — ro|| and we have introduced the Gaunt coeffi- 
cients (see eq. (31) in the text) 

CjK;lm,l'm' = / dn^YjK{k)Y;jk)Y,rn'{k) (A.6) 

Inserting the identity of eq. (A.5) in the partial wave expansion of the greenian 
as given in eq. (A. 3) we arc able to rewrite this expansion in the same form of 
eq. (29) that we have used in the text to derive the quantization conditions, 
i.e. 



9k{r-fo;q)^qJ2jii(lT^) yim{r)ni{qro) ^(^^(fo) 

Im 

+ H ji{qr) Yim{r) rim,i'm'{k,q) ji'iqro) Yi'm'ifo) (A. 7) 

Iml'm' 

where r < Tq < nL ^ 0, and we have introduced the structure coefficients 

^lm,l'm'{k, q) = AtT i^~'^"^i~'^DjK{k, q) CjK-lm,l'm' (A. 8) 

JK 

B Structure coefficients calculation 

In this appendix we derive a mathematical expression useful for the numer- 
ical calculation of the structure coefficients. From eq. (A. 8) we realize that the 
non trivial part of this problem consists in the calculation of the reduced struc- 
ture coefficients. In the following we are going to derive an expression of the 
reduced structure coefficients, different from that already given in eq. (A.4), 
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in order to handle a formula suitable for the numerical evaluation. We start 
again from the greenian definition 



\R-nL\ 



ik-nL 



R-nL 



fB.ll 



having the aim to express it in the reciprocal space. To this end we recall the 
following well known identities 



(2vr) 



2n 



Pm 



neZ3 ^ 



1 



£-»o+ (27r)3 



dp 



1 e*«^ 



jp' — {q^ -\- ie) (47r) x 
that allow us to write the following chain of equalities 



(B.2) 
(B.3) 

(B.4) 



9kA Q) 



1 



(27r)3e^0+ ^ 



lim 



Ak-nL 



neZ3 



Qip-{R—nL) 

p2 _ (•g2 _|_ ^g-J 



2 ^ipR 

lim / dp 



E 



^i(k—p)-nL 



(27r)3e-.0+i ^ p2-^q2 + i^^ ^^^^ 
gip-R I 



^hm J dp ^2 _ ^ (^27r)3 



dx E (^(f-nL)e^M-^ 

neZ3 



1 



Pm 



^ip-R 



1 



= - lim 

£-»0+ 



p^-{q^ + ie) (27r)3 



/ 



— + ie) 



S{k + pm-p) 



1 ^ e 



i{k+pm)-R 



L't:ik+Pmr-q' 
i.e., the greenian expression in the reciprocal space is given by 



(B.5) 



_^ 2 ^iqm-R 

9k{R;Q) = -71E 



- 27r ^ 



(B.6) 



We are now ready to derive an expression for the KKR reduced structure 
coefficients in the reciprocal space. To this end we need to recall the identity 
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expressing a plane wave as an expansion in spherical harmonics 

^n^R ^ ATlY^l^jMruR) YiUqm) Yim{R) (B.7) 



Im 

so that the greenian can be expanded as follows 

In order to reproduce the same structure of eq. (A. 3) we make the following 
observation 

, / V 1 cos(qR) 1 cos(qR) 

, ^ , 1 cos(qR) (7 / / N 

= 9dR-,<l) - + ^cot(gi?) jo{qR) (B.9) 

inserting the last identity in eq. (B.8) we end up with 

9kiR-^ Q) = -^^^^ + E DimiK q; R)ji{qR) YUR) (B.IO) 

Im 

where we have obtained 

Dim[k, q; R) = — — ^ — + ^= cot(gi?) d^o dmo (B.ll) 

This expression requires some comments. First we want to note that 

Ytm{qm) = Yi^[^J q^ = (g^, g^, g^) g'^ = (gj^, -g^, g^) (B.12) 

so that, being the rest of the sum argument a function of the modulus of g^ 
only, we can write 

Di^[k, g; i?) = — — }^ ^ + -y= cot(gi?) (B.13) 

The second observation concerns the i?-dependence of the reduced structure 
coefficients. Indeed, by the comparison of eq. (A. 4) with the previous equa- 
tion we learn that this functional dependence is fictitious. This can be also 
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understood by a careful analysis of eq. (B.IO). We know that the term 

, cos(qR) ,^ 

noiqR) = (B.14) 

satisfies by its own the equation 

(a + q^) noiqR) = -^S{R) (B.15) 

so that the remaining terms of the sum in eq. (B.IO) must be regular solutions 
of the homogeneous Helmholtz's equation. A general solution of this kind can 
be expressed as a linear combination of the form 

ki^-^ Q)=E DimiK q)jiiqR) YimiR) (B.16) 

Im 

where the Dijn{k,q) are constants by respect to R. We also learn from text- 
books on mathematical functions that near the origin the spherical Bessel 
functions deal as 

Ji{qR)^iqRy (B.17) 

The fictitious it!-dependence of the reduced structure coefficients as given in 
eq. (B.13) can be eliminated taking the limit for R that goes to zero and using 
the previous equation. 

One has to consider eq. (B.13) as a properly regulated form of the KKR 
reduced structure coefficients and think to the following 

Di^{k,q) = — Yj^2^— — + lim - B.18) 

qL^t^ Tm-q V47r R^^R 

as a formal expression of the same objects. 
B.l Ewald's sums 



In ref. [13] Kohn and Rostoker considered a method particularly con- 
venient from the numerical point of view to evaluate the reduced structure 
coefficients. They pointed out that, following a prescription due to Ewald [20], 
the sum 

S{x) = E ) ^'-"(^-) (B.19) 

approximate the needed result S{oo) with an exponential vanishing error. In 
ref. [21] the original observation of Kohn and Rostoker was further refined by 
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Ham and Segall. They first recall two identities both due to Ewald; the first 
one is 



J,q R—fiL 



R-nL 



hm ^ / dx e-(^-"^^)^-^+^ 



(B.20) 



The second identity can be proved as done before to express the greenian in 
the reciprocal space starting from its expression in the direct space and states 



^ ^-{R-nLfx'^+ik-inL-R) ^ iV^T ^^ipm-R-^ 
neZ3 



(LxY 



Pm = (B.21) 



We now start again from the expression of the greenian in the direct space as 
given in eq. (B.l) and write, by using the first identity (B.20), 



1 -■ ^ f - jk.fiL-{R-HLfx^+S^ 



lim ^ dx e 



(B.22) 



and 



, where 77 is a 



then we split the integration domain in 0, 2 
positive arbitrary constant. The greenian can be thus re-expressed as the sum 
of two terms 



9dR-^ q) = 9\{R\ q) + ^K^; q) 



(B.23) 



where, performing the integration in the first term and using the identity of 
eq. (B.21) in the second one, we obtain 



9\iR-^q)-jjX 



1 °° 

2 



(B.24) 



These scries arc absolutely convergent for any finite 77 > and each term is an 
analytic function of q throughout the complex plane except at the simple pole 
— Expanding termwise both and in spherical harmonics by respect 

— * 

to i?, taking the limit — > 0, and comparing the result with the definition of 
the reduced structure coefficients, we find 



DUk, q) = DUk, q) + ^^(fc, q) + DL(^, q) 



(B.25) 



where 
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2 

1m 



2 

4x^ 



2 



(B.26) 



In order to make more explicit the dependence of the reduced structure coef- 
ficients upon the volume let us define 



dim{0,Q)=LDim{k,q) 



^dUe,Q) + dUO:Q) + dUO,Q) 



(B.27) 



and rewrite eqs. (B.26) as follows 



{'Kny^'-^'YUn) I dxx 

^ n6Z3-{o} ^ 



1— °° 1 



(B.28) 



where we have substituted 



Q 



Lq 

— * 

e 



Qr, 



m + 



27T 



2^) ^ 



(B.29) 
(B.30) 



B.2 Incomplete Gamma Function 



The computation of the reduced structure coefficients is complicated from 
the numerical point of view by the presence of an integral in the definition of 
d'i^{9,Q). In order to simplify the numerical task let us recall the definition 
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of the incomplete gamma function: 

oo 

r(a,y) ^ j dx a;"-^e-^ (B.31) 
y 

using the incomplete gamma function, d}^{6, Q) can be rewritten in the fol- 
lowing form 



The great advantage of this expression with respect to the one given in eq. (B.28) 
is that the incomplete gamma function can be computed numerically using a 
continued fraction representation: 

r(«, y) = ^ I i-g — (B.33) 

^ — 

X H ij- 

1 _L ^ 



3-a 
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